Theoretical study of molecular electronic excitations and optical transitions of Ceo 
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We report results on ab initio calculations of excited states of the fuUerene molecule by using 
configuration interaction (CI) approach with singly excited determinants (SCI). We have used both 
the experimental geometry and the one optimized by the density functional method and worked 
with basis sets at the cc-pVTZ and aug-cc-pVTZ level. Contrary to the early SCI semiempirical 
calculations, we find that two lowest ^Tiu <— ^Ag electron optical lines are situated at relatively 
high energies of ~5.8 eV (214 nm) and ~6.3 eV (197 nm). These two lines originate from two 
^Tiu ^ ^Ag transitions: from HOMO to (LUMO+1) (6/i„ 3tig) and from (HOMO-1) to LUMO 
{Whg — > 7tiu)- The lowest molecular excitation, which is the l^T2g level, is found at ~2.5 eV. 
Inclusion of doubly excited determinants (SDCI) leads only to minor corrections to this picture. We 
discuss possible assignment of absorption bands at energies smaller than 5.8 eV (or A larger than 
214 nm). 

PACS numbers: 31.25.-v, 36.20.Kd, 81.05.Tp 



I. INTRODUCTION 

The buckminsterfuUerene molecule Ceo attracts much 
attention of theoreticians since its discovery in 1985 
Indeed, the role of the fuUerene molecule in quantum 
chemistry is unique because of its highest group (Ih) of 
symmetry leading to a high degeneracy of molecular or- 
bitals (MOs). The MOs of Cgo are classified according to 
irreducible representations (irreps) of 7^, which are a, ti, 
t2, g, and h. (The dimensions of the irreps are 1, 3, 3, 4, 
and 5, correspondingly.) For example, the highest occu- 
pied molecular orbital (HOMO) is a fivefold electron shell 
of hu symmetry while the lowest unoccupied molecular 
orbital (LUMO) is a threefold level of tiu symmetry. 

High degeneracies of molecular orbitals imply a com- 
plex electronic structure [1, 0, H, @]. However, since all 
electron shells of the neutral Ceo molecule are completely 
filled, the ground state is of ^Ag symmetry and the many- 
electron states of the molecule are not revealed. The sit- 
uation changes when the fuUerene molecule is excited. 
Then a few electrons of Ceo are promoted to a higher 
energy, and some of electron shells become open. The 
electrons of the open shells exhibit correlated behavior 
and form many electron levels called molecular terms 
[3, 0, H, H, , which are also classified according to the 
irreducible representations of the icosahedral symmetry 
{A, Ti, T2, G, and H). Thus, a careful study of these cor- 
relations is needed for calculations of optical transitions 
[ill ll^ and molecular excitations. It is worth mention- 
ing that intramolecular correlations play an important 
role in alkali doped fuUerites, leading to superconduc- 
tivity, magnetic behavior and metal-insulator transitions 



TABLE 1: Selected molecular excitations calculated with 
semiempirical CI calculations and in the present work. Ener- 
gies are in eV. Ncsf stands for the number of configuration 
state functions used for CI. 



Ref. 


[9] 


n 


flO] 


this work 


method 


PPP 


QCFF/tt 


CNDO/S 


ah initio SCI 


Ncsf 


134 


266 


900 


2857 




planar 3D 










2.456 2.232 


2.06 




2.549 




3.603 3.381 


3.21 




4.700 




3.635 3.391 


3.15 




4.771 


i'ri„ 


4.227 4.000 


4.08 


3.40 


5.796 


2 Tin 


4.750 4.681 


4.53 


4.06 


6.335 
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Baa. 

Computations of many electron states of Ceo is a chal- 
lenge for the modern quantum theory (isj . In the lit- 
erature there are few works devoted to the problem 
of calculations of molecular excitations of the fuUerene 
0, B S 113], but to the best of our knowledge all of 
them are semiempirical and are based only on singly ex- 
cited configuration interaction (CI). We reproduce some 
of their computed levels in TablelH The calculation of Ne- 
gri et ai, Ref. employed a semiempirical (QCFF/tt) 
method, and took into account 196 configuration state 
functions (CSF). The first excited state {l^T2g) was 
found at 2.06 eV, while l^Tiu and 2^Tiu levels at 4.08 and 
4.53 eV. A PPP-CI (Pariser-Parr-Pople Configuration In- 
teraction) calculation which includes 134 CSF finds first 
two ^Tiu levels at 4.00 eV, 4.68 eV and 6.69 eV, with 
the first excited ^Tzg level at 2.23 eV. Finally, CNDO/S 
calculations of Braga et al. based on 808 and 900 configu- 
rations predicts that l^Tiu is lowered to 3.4 eV, while the 
three most intense transitions are at 4.38, 5.24 and 5.78 
eV. Unfortunately, all of the CI calculations reported in 
Refs. 0, BBIi^l are crucially dependent on a number of 
approximations for matrix elements and therefore cannot 
be considered as a true ah initio approach. 

In the present study in comparison with Refs. 0, H, 
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[lO| we employ the ab initio treatment, which represents 
a new level of CI calculations of the fuUerene molecule. 



II. METHOD OF CALCULATION 
A. Configuration Interaction 

A well known recipe for correlations is the method of 
configuration interaction (CI) The basic idea is to 

diagonalize the A^-electron Hamiltonian in a basis of A^- 
electron functions (Slater determinants). In principle, 
the full CI method with infinite number of the molecular 
orbitals provides an exact solution of the many-electron 
problem. In practice, even for small molecules and mod- 
erately sized one-electron basis sets, the number of A^- 
electron determinants is enormous. To avoid it, we first 
introduce an active space which comprises a set of highest 
occupied molecular orbitals (HOMOs) and a set of low- 
est unoccupied molecular orbitals (LUMOs). The space 
of the MOs is divided into three subspaces: the inactive, 
the active and the external orbitals. The inactive orbitals 
are double occupied, while the external orbitals are unoc- 
cupied. Within the chosen active space one can consider 
various approximations to the complete active space CI 
matrix by truncating the many-electron trial function at 
some excitation level. For example, singly excited CI 
(SCI) approach has been proven to be adequate as a first 
approximation for the neutral molecule because all elec- 
tron shells of Ceo are filled up. In the following for large 
active spaces (AS) we use SCI, but also consider singly 
and doubly excited determinants (SDCI) for smaller AS, 
and a complete active space (CAS) CI method fl^ for 
important active spaces of two molecular shells (8 MOs). 
The computer code for the CI calculations is the modified 
and extended version of the earlier program, Ref. 



B. Computational details 

The CI calculations have been carried out with the set 
of molecular orbitals (MOs) obtained from the restricted 
Hartree-Fock (RHF) self-consistent-field calculation [l^ 
of the neutral fuUerene molecule. In the following we have 
used two sets of coordinates of the Cgo molecule. First 
set was obtained in Ref. fl7l (DFT/B3LYP method). 
Two different C-C bond lengths were i?5 = 1.4507 A 
(C-C bond in pentagons) and Rq = 1.3906 A (short C- 
C bond in hexagons). The^second set corresponds to 
the experimental geometry jlS] with Rc, — 1.448 A and 
Re = 1.404 A. 

We have used our original RHF computer program [ll| , 
which uses the resolution of the identity ^3] (M) method 
for calculation of the two-electron Coulomb integrals. 
The series of the RI (auxiliary) basis sets, HCNO.x, were 
designed for the Rl-convergent calculations the RHF to- 
tal energies and the electronic properties of molecules 



TABLE II: Electron shells of Ceo- Emo is the one-electron 
Hartree-Fock MO energy (in eV). The basis set is aug-cc- 
pVTZ(-f). Etot = -61832.163 a.u. 





MOs 


symmetry 


Emo 


LUMO-HIS 


244-248 


Shn 


4.303 


LUM0-H7 


239-243 


Uh, 


3.855 


LUM0-fl6 


235-238 




3.794 


LUMO-HIS 


232-234 


10*2^, 


3.300 


LUMO-HM 


229-231 




3.280 


LUMO-M3 


228 


Gag 


3.161 


LUMO-hl2 


225-227 




3.115 


LUMO-hll 


220-224 




2.557 


LUMO-HIO 


215-219 


13/lg 


2.418 


LUMO-h9 


211-214 


759 


2.254 


LUMO-H8 


206-210 


12hg 


2.121 


LUMO-f7 


203-205 


9t2u 


1.994 


LUMO-F6 


200-202 




1.774 


LUMO-h5 


196-199 


TQu 


1.751 


LUMO-H4 


195 


Scig 


1.496 


LUMO-f3 


190-194 


llh. 


1.398 


LUMO-H2 


187-189 


lt2u 


1.150 


LUMO-Hl 


184-186 


3ii9 


0.947 


LUMO 


181-183 


7tiu 


-0.817 


HOMO 


176-180 




-7.834 


HOMO-1 


171-175 


lO/lg 


-9.627 


HOMO-2 


167-170 


659 


-9.890 


HOMO-3 


163-166 




-12.589 


HOMO-4 


160-162 


&t2u 


-13.050 



containing hydrogen, carbon, nitrogen and oxygen ele- 
ments as discussed in [l^. Here x indicates (but not 
exactly equal to) the accuracy (in atomic units) of the 
total RHF energy: smaller x implies more RI functions 
and more accurate RI basis set. The accuracy of the RI 
basis set is almost independent of the choice of molec- 
ular basis set. Having performed test calculations with 
Pople's 6-31G* MO basis set ,21, ,22] and the HCNO.x RI 
series we have found out that a reasonable convergence 
is achieved for the HCNO.OOl RI basis set, which is used 
for all further calculations. 

Throughout the paper we have adopted short nota- 
tion aug-cc-pVTZ(-f) for Dunning's augmented correla- 
tion consistent polarizable valence triple zeta molecu- 
lar basis set aug-cc-pVTZ [2^ where all polarization 
functions of /-symmetry are excluded from the set. 

III. RESULTS AND DISCUSSION 

A. SCI excitation spectrum of Ceo 

The main results are shown in Tables [Til and IIIIl The 
(HOMO-4) to (LUM0-H18) energies of Ceo calculated 
with the aug-cc-pVTZ(-f) basis set are quoted in Table [III 
We then define a broad active space spanned by these 
MOs and occupied with 42 correlated electrons and use 
it in our SCI calculations. Table IIIIl displays the result- 
ing spectrum of lowest molecular excitations. In order to 
evaluate the influence of two bond lengths of the fuUerene 



3 



TABLE III: Lowest excitations of Ceo- OG refers to the 
bond lengths of Ceo optimized by DPT, Ref. EG to the 
experimental ones, Ref. fisl . 





'Tig 


'T2, 




^7^1. 


^Gg 


^Gg 


'Gu 


OG 2.549 3.041 


3.229 


3.287 3.416 


3.417 


3.538 


3.595 


EG 2.384 


2.877 3.051 


3.100 


3.239 


3.235 


3.357 


3.464 


1 lu 






'Hu 




Gu 


'^Hu 


-t2ti 


OG 3.917 3.996 4.079 4.304 4.611 4.700 4.771 4.960 


EG 3.771 


3.801 


3.906 


4.154 


4.447 4.534 4.588 


4.822 




Gu 


"^Gg 


' Hu 


^Hu 


' Hu 


^Tiu 


J lu 


OG 5.163 


5.291 


5.355 


5.454 


5.598 


5.732 


5.796 


5.972 


EG 5.023 


5.157 


5.300 


5.295 


5.445 


5.662 


5.667 


5.849 


'Hu 




"^Gu 




^ Hu 


Gu 


'Hu 


^Hu 


OG 5.986 


6.001 


6.019 


6.124 


6.141 


6.208 


6.236 


6.264 


EG 5.823 


5.872 


5.873 


6.041 


5.954 


6.021 


6.188 


6.218 


"'"Tin 


'Gu 


J lu 








' Hu 




OG 6.335 6.337 6.340 6.412 6.422 


6.444 


6.501 


6.518 


EG 6.182 


6.216 


6.206 


6.340 


6.275 


6.360 


6.395 


6.400 



molecule, we have performed calculations with two dif- 
ferent geometries of Ceo: optimized by DFT [17] and the 
experimental one [l^. Inspection of Table Hill shows that 
the order of levels is mainly conserved, while a typical 
numerical deviation is ~ 0.2 eV. From this we conclude 
that a possible small change of the Ceo bond lengths do 
not change appreciably its excitation spectrum. 

The most striking feature of the calculation is the rela- 
tively high energy position of the two lowest ^Ti„ excita- 
tions (given already in Table [T| with oscillator strengths 
1.02 and 0.16, respectively. For a review on oscillator 
strengths in semiempirical calculations see Ref. [l^]. In 
the following we limit ourselves to the optimized geome- 
try [17] and consider mainly the problem of convergency 
and adequacy of the SCI calculation leaving discussion 
and conclusions until Sec. IIVI 

B. Basis set dependence 

To study basis set dependence of the SCI excitation 
spectrum we narrowed active space to 167-219 MOs, see 
Table [III This active space consists of 14 electron shells 
(from HOMO-2 to LUMO+10) filled by 28 correlated 
electrons. Such calculations require less computer time in 
comparison with the calculation with the active space of 
Table nil reported in Table HIH but lead to a less accurate 
molecular excitation spectrum. 

The results of such SCI calculations are presented in 
Table IIVI Notice that the order of the levels is almost 
the same for the group of 6-31G*, cc-pVDZ, cc-pVTZ 
and cc-pVTZ(-f ) basis sets. The deviations of the molec- 
ular excitation energies are rather small, of the order of 
0.1 eV. The molecular energies are close to a convergence 
and the basis f-functions have a small effect on calcu- 
lated values. On the other hand the inclusion of diffuse 
functions of aug-cc-pVDZ(-f) and aug-cc-pVTZ(-f) leads 
to a rearrangement of electronic levels, Table IIVI which 
indicates that aug-cc-pVTZ(-f) is the best set for SCI 



TABLE IV: Energy spectrum of Ceo with various basis sets. 
Basis sets are cc-pVDZ [DZ], aug-cc-pVDZ [aDZ)], cc-pVTZ 
[TZ)], cc-pVTZ without the /-functions [TZ(-f)], and aug- 
cc-pVTZ(-f) [aTZ(-f)], see text for details. 



r 


6-31G* 


DZ 


TZ 


TZ(-f) 


r 


aDZ 


r 


aTZ(-f) 




2.593 


2.575 


2.582 


2.573 


Tig 


2.900 


T2g 


2.669 


Tig 


3.114 


3.071 


3.053 


3.048 


Tig 


3.261 


3r7-i 
Tig 


3.107 


1 m 

T2a 


3.322 


3.276 


3.249 


3.248 


1 m 

T2g 


3.325 


1 m 

Tig 


3.254 


3 rr 

H, 


3.384 


3.336 


3.312 


3.308 


3 rr 
Hg 


3.474 


3 rr 
Hg 


3.390 


1 m 


3.522 


3.467 


3.434 


3.432 


Tig 


3.541 


1 r-r-i 

Tig 


3.439 




3.540 


3.481 


3.447 


3.445 


3 

Lrg 


3.547 


Lfg 


3.477 


Gg 


3.654 


3.594 


3.558 


3.556 


Lrg 


3.649 


Lrg 


3.578 


Gu 


3.672 


3.675 


3.700 


3.685 


lu 


4.154 


3 

Gu 


3.995 


3rT-i 

-t lu 


3.899 


3.884 


3.884 


3.875 


1 2u 


4.196 


3rj-i 

lu 


4.069 


J-2u 


4.002 


3.993 


3.995 


3.987 


Hg 


4.288 


1 2u 


4.132 


1 TT 

Hg 


4.231 


4.159 


4.112 


4.111 


3 

Gu 


4.355 


1 

Hg 


4.136 


3 rr 


4.383 


4.365 


4.365 


4.356 


3 rr 
tLu 


4.754 


3 rr 
tLu 


4.541 


1 rri 


4.650 


4.632 


4.625 


4.618 


1 2u 


4.889 


1 2u 


4.751 


1 y~i 


4.831 


4.802 


4.779 


4.777 


1 

Gu 


4.897 


1 

Gu 


4.849 


1 rr 


4.862 


4.835 


4.816 


4.811 


1 rr 
tlu 


4.976 


1 rr 

tiu 


4.899 


J-2u 


5.120 


5.075 


5.053 


5.045 


1 2u 


5.500 


3rj-\ 

± 2u 


5.241 


3 

Git 


5.402 


5.337 


5.293 


5.289 


3 

Gu 


5.651 


3 

Gu 


5.350 


3 /"f 


5.443 


5.398 


5.392 


5.373 


1 

Gu 


5.771 


1 

Gu 


5.451 


1 /~f 

Git 


5.514 


5.448 


5.404 


5.400 


3 rr 

-Ou 


5.603 


3 rr 
t±u 


5.529 


3 ZLT 


5.576 


5.524 


5.490 


5.486 


1 rr 
tLu 


5.722 


1 rr 

tiu 


5.652 


1 TJ 


5.764 


5.704 


5.661 


5.659 


Irji 

lu 


6.116 


i lu 


5.880 


1 lu 


5.873 


5.826 


5.793 


5.785 


3 rr 
tLu 


6.221 


3 

Lfg 


5.919 


3 rr 
J^u 


5.940 


5.859 


5.810 


5.797 


1 2u 


6.236 


i 2u 


6.071 


3/7^ 
-t lu 


6.00( 


0.992 


0.94D 


0.93/ 


Ug 


0.249 


3 rr 

tiu 


O.Ulb 




6.087 


6.032 


5.993 


5.987 


3t^ 
lu 


6.280 


3t^ 
lu 


6.092 


Gu 


6.116 


6.067 


6.030 


6.022 


^Hu 


6.254 


'Hg 


6.098 


^T2u 


6.165 


6.096 


6.048 


6.046 


Gu 


6.281 


Gu 


6.103 


''T2g 


6.194 


6.113 


6.054 


6.042 


''T2g 


6.543 


^Hu 


6.198 


Hu 


6.299 


6.246 


6.202 


6.199 


Tiu 


6.577 


Gu 


6.264 


Gu 


6.363 


6.314 


6.274 


6.271 


Gu 


6.954 


''T2g 


6.325 


Gu 


6.445 


6.386 


6.346 


6.340 


''Hg 


6.960 


Tiu 


6.371 


-[ lu 


6.457 


6.402 


6.360 


6.356 


J lu 


7.118 


J lu 


6.413 


Tiu 


6.462 


6.407 


6.365 


6.362 


Tiu 


7.122 


Gu 


6.452 



calculations. 



C. Active space dependence 

Here we study the dependence of SCI excitation ener- 
gies on the choice of active space. In Table |V] we keep the 
lower boundary active shell, which is the IIOMO-4 {6t2u) 
level unchanged and systematically increase the upper 
electron boundary. Inspection of Table |V] shows that we 
are close to a convergence for the broad active space of 
160-248 MOs, shown in Table [U Therefore, we do not 
expect drastic changes if the active space is increased 
even further. Notice however, that the optically active 
l^Tiu and 2^Ti„ levels monotonically decrease with the 
increase of active space. 
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TABLE V: Dependence of selected energy levels on cho- 
sen active space. The active space consists of the MOs be- 
tween lower and upper MOs. The lower MO is 160, which 
corresponds to the HOMO-4 (6*21,) level. The upper MO 
changes from LUMO-f 1 to LUMO-l-18. The basis set is aug- 
cc-pVTZ(-/). 

upper MO 186 189 194 199 202 205 2T0 2T4~ 
l^Tag 2.889 2.889 2.889 2.889 2.887 2.833 2.713 2.712 
l^Gu 4.808 4.808 4.808 4.808 4.807 4.801 4.758 4.758 
l^Hu 4.898 4.898 4.898 4.898 4.898 4.891 4.836 4.836 
l^Tiu 6.230 6.230 6.230 6.230 6.230 6.212 6.028 6.026 
2^Tiu 7.210 7.210 7.194 7.172 7.112 6.936 6.421 6.419 

upper MO 219 224 227 231 234 238 243 248~ 
l^Tag 2.633 2.633 2.623 2.618 2.600 2.599 2.584 2.549 
l^Gu 4.730 4.730 4.727 4.721 4.719 4.719 4.715 4.700 
l^Hu 4.802 4.802 4.798 4.793 4.791 4.790 4.785 4.771 
l^Tiu 5.862 5.862 5.862 5.860 5.839 5.838 5.816 5.796 
2^Tiu 6.354 6.354 6.352 6.352 6.348 6.348 6.343 6.335 



D. Beyond SCI approximation for small AS 

We recall that all the results reported earlier have 
been obtained under approximation that only singly ex- 
cited configurations (SCI) are taken into account. It is 
therefore instructive to estimate the accuracy of this as- 
sumption. Since the active space shown in Table |TT] is 
very big, we study this problem by considering fewer 
active MOs. We limit ourselves to four most im- 
portant active spaces: HOMO-I-LUMO {6hu + 7ti„), 
HOMO-f (LUMO-fl) (6/i„ -I- 3tig), (HOMO~l)-hLUMO 
{10hg + 7tiu), and (HOMO~l)-H(LUMO-f l)(10/ig-f 3iig). 
All of them consists of 8 active MOs (five from h and 
three from ti), which accommodate 10 correlated elec- 
trons. For each of the active spaces we have performed 
four calculations: SCI, singly and double excited CI (SD- 
CI), singly, doubly and triply excited CI (SDT-CI), and 
with all possible excitations. The latter is so called com- 
plete active space (CAS-CI) calculation [l^. In Tables 
I VIIIIXI we give the results for the four cases, respectively. 
Notice that the CAS-CI results are very close to that of 
SDT-CI, while both SCI and SDCI molecular energies 
miss the precise values by a typical error of ~ 0.1 eV. 
Furthermore, as a rule SDCI energies overestimate the 
real values. It is very probable that the same conclusions 
are applicable to a calculation with a large AS. Indeed, 
our SDCI calculations [25] with a large AS of 186-248 
MOs unambiguously indicate an increase of ^Ti„ ener- 
gies in comparison with the SCI values: the energy of 
the first transition changes from 5.997 to 6.291 eV, while 
the energy of the second from 6.395 to 7.154 eV. Thus 
the SDCI results underline the importance of SDT-CI 
calculation for the Ceo molecule. 

It is worth noting that two lowest ^Ti„ ^ ^Ag transi- 
tions can be found already in Tables IVlTl and IVIIII Their 
energies are quite high, 6.44 eV (6/i„ 3tig) and 6.96 eV 
(10/ig — » 7tiu)- If both of the ^Tig levels are included in 
a CI calculation, they interact with each other, because 



TABLE VI: Lowest excitations of Ceo, calculated with 10 
correlated electrons in the 6/iu -I- 7iiu limited active space 
(see text for details). 



excitation level 


s 


SD 


SDT all (CAS) 


1 /I 














^T2g 


3.312 


3.383 


3.338 


3.338 


'Tig 


3.482 


3.551 


3.511 


3.511 


'T,g 


3.590 


3.661 


3.622 


3.623 


'Hg 


3.640 


3.699 


3.651 


3.651 


]Tlg 


3.691 


3.761 


3.713 


3.713 


Gg 


3.709 


3.768 


3.726 


3.726 


'Gg 


3.769 


3.833 


3.791 


3.791 


'Hg 


4.298 


4.337 


4.281 


4.280 



TABLE VII: Lowest excitations of Ceo, calculated with 10 
correlated electrons in the 6/iu -I- 3tig limited active space (see 
text for details). 



excitation level 


S 


SD 


SDT all (CAS) 


M 














J-2u 


5.188 


5.322 


5.239 


5.239 


-[ lu 


5.234 


5.368 


5.297 


5.298 


Gu 


5.410 


5.544 


5.470 


5.471 




5.461 


5.595 


5.529 


5.530 




5.550 


5.684 


5.616 


5.617 


Gu 


5.594 


5.728 


5.661 


5.662 




5.692 


5.826 


5.760 


5.761 




6.443 


6.577 


6.438 


6.439 



TABLE VIII: Lowest excitations of Ceo, calculated with 10 
correlated electrons in the lOhg + 7tiu limited active space 
(see text for details). 



excitation level 


S 


SD 


SDT all (CAS) 


^9 














Gu 


5.218 


5.497 


5.331 


5.334 


-[ lu 


5.233 


5.512 


5.371 


5.374 




5.440 


5.719 


5.562 


5.565 


'^Hu 


5.457 


5.737 


5.596 


5.600 


^T2u 


5.548 


5.828 


5.679 


5.682 


Gu 


5.580 


5.860 


5.716 


5.720 


^Hu 


5.605 


5.884 


5.737 


5.741 




6.997 


7.277 


6.959 


6.960 



TABLE IX: Lowest excitations of Ceo, calculated with 10 
correlated electrons in the lO^g -I- 3tig limited active space 
(see text for details). 



excitation level 


s 


SD 


SDT all (CAS) 


Mg 














'Hg 


7.072 


7.110 


7.075 


7.075 


'T2g 


7.116 


7.154 


7.125 


7.125 


'Gg 


7.174 


7.224 


7.194 


7.194 


'Tig 


7.192 


7.236 


7.209 


7.209 


'T2g 


7.299 


7.348 


7.320 


7.320 


'Tig 


7.337 


7.383 


7.353 


7.353 


'Gg 


7.388 


7.442 


7.415 


7.415 


'Hg 


7.868 


7.898 


7.859 


7.859 
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FIG. 1: Comparison between the experiment data on optical 
density, Ref . jllJ, and the calculated spin singlet lines of Ceo- 
The dashed line is the contribution from two Tiu transitions 
broadened with the line width of 0.5 eV. The MO basis set is 
aug-cc-pVTZ(-f), the bond lengths are taken from Ref. [TtI. 



they are of the same symmetry. Nevertheless, the lowest 
l^Tiu level stays predominantly of the 6hu — > 3iig origin, 
while the 2^Tiu level is mainly of the lOhg Itiu charac- 
ter. Notice that the lowest part of molecular excitations 
of Ceo stems from the HOMO-^LUMO ['oK + ltiu) active 
space, Table Ell 



IV. DISCUSSION AND CONCLUSIONS 

We have found that the lowest electron optical transi- 
tions are at ~ 5.8 eV (214 nm) and ~ 6.3 eV (197 nm), 
which are at disagreement with the results obtained with 
the early semiempirical CI calculations [1, [l3|. We 
have analyzed the influence on the SCI excitation spec- 
trum from the approximations used in our calculation 
and concluded that all of them are unlikely to change 
the computed values substantially. Thus, according to 
our calculations, the first allowed optical transition of 
Ceo lies very close to the characteristic interstellar ab- 
sorption at 217 nm 

From the ab initio study we come to the conclusion 
that early semiempirical configuration interactions of the 
fuUerene molecule, Refs. result in systematically 

smaller excitation energies, see Table ID This was also 
the case for the molecular ion Cgo , for which the ^Ag 
spin singlet ground state was obtained [27j . As shown in 
Refs. [5, 6] the ground state is rather the ^Tig spin triplet 
in accordance with Hund's rules. 

Calculated oscillator strengths of two lowest optical 
transitions to l^Ti„ and 2^Ti„ states in the largest AS 
(Table HH are 1.02 and 0.16, respectively. QCFF/tt 
method employed by Negri et al. Q resulted in 0.61 



and 0.41 for oscillator strengths of two lowest transitions 
while Braga et al. jlOi] [CNDO/S method] obtained 0.08 
and 0.41. As we remarked before our ab initio energies 
of two lowest ^Tiu transitions are very different from the 
semiemperical ones and a direct comparison of the oscilla- 
tor strengths is probably premature. On the other hand, 
early semiemperical calculations give only a vague cor- 
respondence with the experimental optical density pT| . 
see Fig. 14 of Ref. [l^l . In order to achieve a good com- 
parison Westin et al. used a screened response oscillator 
strength distribution calculated with an adjustable pa- 
rameter V [23|. Under such conditions the theoretical 
curve reproduces the three experimental bands shown in 
Fig. 1 . Our present ab initio study accounts only for the 
6 eV band. Fig. 1. Thus, two other bands are probably of 
the electron- vibronic origin [28,] . The situation is not un- 
common for complex molecules and it occurs for example 
at the ^B2u *~ ^^ig (3200 A) transition of naphthalene 
as discussed in Chapter II7 of Ref. jl^ . 

The optical spectrum of Ceo also shows some similari- 
ties with that of benzene [2^, [23 . Lowest excited states 
of benzene are spin triplets, while the first electronically 
allowed optical transition (l^i?i„ ^ ^jdls) is situated 
at relatively high energy of ~7 eV [2^ |30| . In benzene 
two lowest absorption bands at 4.90 and 6.20 eV were 
unambiguously assigned to the excited states of l^B2u 
and l^Biu, respectively [1^, [13]. These bands are elec- 
tronically forbidden and occur only due to the electron- 
vibronic coupling, which underlines the importance of 
the Herzberg- Teller mechanism ^] for the molecule. It 
is plausible that the same scenario applies to the Ceo 
fuUercne, and optical lines at smaller energies are caused 
by vibronic coupling to numerous electron states in these 
energy regions. Fig. 1. Indeed, there are indications that 
optical bands observed in the low region of the absorp- 
tion spectrum of Ceo can be interpreted as false origins of 
states of the Tig symmetry [3l|. It is worth noting that 
singlet and triplet states of Ceo can be distinguished by 
low energy electron spectroscopy [32j . 

Unfortunately, the electronic spectrum of Ceo is not so 
well understood as benzene's, which has a long history of 
calculations and interpretations [1^, [2^ . Although in all 
our SCI and SD-CI calculations the l^Ti„ level remains 
at high energy of '^6 eV, we think that there is still a pos- 
sibility for it to be lowered if a calculation with a higher 
level of excitations within the CI approach is performed 
(at least SDT-CI). Interestingly, this is the case for the 
first allowed optical transition (l^i?i„ ^ ^^ig) in ben- 
zene. The CASSCF calculation of CeHe (6 tt electrons 
distributed among 12 tt MOs) places the l^Eiu level at 
a high energy of 8.77 eV, while a multireference M0ller- 
Plesset perturbation theory (MRMP) on the CASSCF 
states lowers this energy by 1.84 eV to 6.93 eV (Table 2 
of Ref. [2^ ) . The last value is in perfect agreement with 
the experiment (6.94 eV, Ref. JiOj). Notice that MRMP 
is known for effectively accounting for high excitations. 

In conclusion, our S-CI and SD-CI calculations in- 
dicate that the first electronically allowed transition 
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<— ^Aig is located at a relatively high energy of 5.8- 
6.0 eV. Our finding opens up a very interesting question 
on assignment of two lowest absorption bands at 3.8 and 
4.8 eV found experimentally [lll,[l4|- There are two ways 
to reconcile our calculations with the experiment. First, 
the l^Tiu excited state can be lowered at the level of 



SDT-CI or more refined CI calculations. The second sce- 
nario is the electron-vibronic coupling (Herzberg- Teller 
mechanism) [2^, which is the case for benzene. Further 
investigations and ab initio calculations are needed to 
clarify this issue. 
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